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This  project  focuses  on  an  integrated  expcrimcntal/modeling  study  of  hydrocarbon  oxidation 
mechanisms  under  conditions  of  extreme  thermal  non-equilibrium.  Large,  and  highly  non- 
equilibrium,  initial  pools  of  important  hydrocarbon  combustion  intermediate  and  electronically 
excited  air  species  are  created  in  a  short  pulsed  fuel-air  plasmas,  created  using  -10-20  nsec 
duration  -  high  (-20  kV)  voltage  pulsers,  capable  of  operation  at  repetition  rates  as  high  as  40-50 
kHz.  The  time  evolution  of  critical  species  and  temperature,  after  application  of  a  single 
nanosecond  discharge  pulse,  or  a  rapid  “burst”  of  pulses,  depending  upon  the  experiment,  arc 
experimentally  determined  using  advanced  laser-based  optical  diagnostic  methods.  Results  arc 
compared  to  theoretical  predictions  from  discharge  and  plasma  chemical  kinetic  codes, 
developed  as  part  of  this  program,  in  order  to  provide  model  validation  and  to  improve 
understanding  of  fundamental  nonequilibrium  plasma  kinetics  in  air  and  air/fucl  mixtures. 

Summary  of  Accomplishments 

The  following  five  major  accomplishments  were  the  direct  result  of  funding  provided  by  this 
program. 

i.  Two  Photon  Laser  Induced  Fluorescence  (TALIF)  studies  of  atomic  oxygen  creation  and 

loss  kinetics  in  air,  methane/air,  and  ethylene/air  nanosecond  pulsed  discharges. 

ii.  Pure  Rotational  Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS)  studies  of  oxidation 

and  heat  release  in  ethylene-air  nanosecond  pulsed  discharges. 

iii.  Single  photon  Laser  Induced  Fluorescence  (LIF)  studies  of  NO  creation  and  lost  kinetics 

in  air,  methane/air,  and  cthylene/air  nanosecond  pulsed  discharges. 

iv.  Pure  rotational  Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS)  studies  of  oxidation, 

heat  release,  and  ignition  in  hydrogen-air  nanosecond  pulsed  discharges. 

v.  Development  of  a  new,  non-empirical  model  for  nanosecond  pulsed  discharge  plasma 

dynamics  and  energy  coupling. 

The  most  important  findings  from  each  of  the  above  will  be  briefly  summarized  below.  More 
detail  can  be  found  in  the  cited  publications. 

Experimental  Apparatus 

Figure  1  shows  a  simplified  schematic  of  a  1-D  pulsed  discharge,  optical  access  “flow”  reactor 
that  was  used  for  the  measurements  performed  for  this  program.  The  channel  is  220  mm  long  x 
22  mm  width  x  10  mm  height,  with  1.75  mm  thick  walls,  and  is  constructed  from  a  single  piece 
of  UV  grade  quartz.  Measurements  have  been  conducted  in  dry  air  and  in  mixtures  of  dry  air 
with  methane,  ethylene,  and  hydrogen.  Flow  rates  of  air  and  fuel  are  controlled  by  mass  flow 
controllers,  at  pressures  in  the  range  40  -  100  Torr  and  flow  velocity  of  approximately  u=0.8 
m/sce,  corresponding  to  flow  residence  time  in  the  discharge  region  of  approximately  75  msec. 
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Two  rectangular  copper  plate  electrodes,  14 
mm  wide  by  65  mm  long,  and  are  rounded  at 
the  comers  to  reduce  the  electric  field 
nonuni  formity,  arc  attached  to  the  outside  of 
the  quartz  channel.  Pulsed  plasmas  are 
produced  by  application  of  ~20  kV  peak 
voltage  -  -10-25  nanosecond  pulses,  using  a 
variety  of  high  voltage  pulsed  power  supplies 
available  in  our  laboratory,  capable  of  operation 
at  maximum  repetition  rates  in  the  range  40  to 
100  kHz.  Laser  beams  are  incident  from  either 
right  or  left,  depending  upon  the  experiment.  A 
variety  of  detectors,  filters,  grating 
spectrometers  /  Optical  Multi-Channel 
Analyzers,  1CCD  cameras  are  available  for 
detection  of  L1F,  TAL1F,  and  CARS  signals,  as 
well  as  for  imaging  discharge  structure. 

Summary  of  Individual  Studies 

i.  Two  Photon  Laser  Induced  Fluorescence  (TAL1F)  studies  of  atomic  oxygen  creation  and 
loss  kinetics  in  air,  methane/air,  and  cthylene/air  nanosecond  pulsed  discharges. 

This  work  is  described  in  detail  in  Publications  2.  9  and  10.  Summarizing  the  experiment  very 
briefly,  the  TALIF  diagnostic  employs  the  well  known  atomic  oxygen  two  photon  allowed  2P  P 
-  3P  ’P  absorption  transition  of  at  -226  nm.  Relative  atomic  oxygen  concentration  data  is  put  on 
an  absolute  scale  using  the  xenon-based  optical  calibration  scheme  ofNiemi,  et  al  [1,2], 

Figures  2  and  3,  which  illustrate  the  most  significant  results,  show  absolute  atomic  oxygen  mole 
fractions  in  air,  (I>  =  1.0  methane/air  (Fig.  2)  and  d>  =  0.5  ethylene-air,  at  60  torr  and 
approximately  300  K  as  a  function  of  time  after  initiation  of  a  single,  20  nsee  discharge  pulse, 
along  with  plasma  chemical  modeling  predictions.  Note  that  the  air  data  in  Figs.  2  and  3  is 
identical,  but  plotted  on  linear  (Fig.  2)  and  log  (Fig.  3)  scales  to  provide  better  comparison  with 
the  fuel  data.  It  can  be  seen  that  in  air  oxygen  atoms  peak  on  a  time  scale  of -10-100  psec  after 
discharge  initiation,  due  to  collisions  with  nitrogen  metastable  molecules  N2(A3X)  +  CL  — ► 
N2(X'E)  +  0  +  0,  which  accounts  for  approximately  75%  of  the  total  molecular  oxygen 
dissociation  produced  by  the  pulsed  discharge.  Atomie  oxygen  loss  in  air  is  due  to  slow  three 
body  ozone  formation,  O  +  O2  +  M  — »  O3  +  M,  which  occurs  on  a  time  scale  of  several  msecs. 

From  Fig.  2  it  can  be  seen  that  in  stoichiometric  methane/air  the  rate  of  atomie  production  is 
essentially  identical  to  that  in  air,  and  the  loss  rate  is  increased  by  a  factor  of  approximately 
three.  As  discussed  in  _,  in  methane/air  at  room  temperature  the  principal  loss  mechanism  for 
atomie  oxygen  is  reaction  with  methyl  (CH3)  radicals,  produced  by  both  electron  impact  and 
dissociating  collisions  of  methane  with  metastable  N2.  Reaction  with  methane,  at  room 
temperature,  is  extremely  slow.  However,  in  ethylene/air  mixtures,  as  shown  in  Fig.  3,  the  rate  of 
decay  of  atomie  oxygen  is  -three  orders  of  magnitude  faster  due  to  processes  such  as  O  +  C,H 

— >  CHj  +  HCO  and  O  +  C  H4  — *  H+  CH^CHO  which  are  fairly  rapid  (rate  coefficient  of  - 


Figure  1.  Schematic  diagram  of  basic  plasma 
assisted  combustion  test  section. 
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k=4.910  cm  /s)  at  -300  K.  Modeling  results  predict  that  rapid  reaction  of  O  atoms  and 
ethylene  initiates  a  series  of  exothermic  low  temperature  oxidation  chemistry,  resulting  in 
substantial  heating  (-500-700  °C)  and,  in  some  eases,  ignition. 


O  atom  mole  fraction 


O  atom  mole  fraction 


Figure  2.  Oxygen  atom  mole  fraction  vs.  time 
after  a  single  high-voltage  pulse  in  air  and  in  a 
methane-air  mixture  at  P=60  torr  and  <1>=1 .0. 


Figure  3.  Oxygen  atom  mole  fraction  vs.  time 
after  a  single  high-voltage  pulse  in  air  and  in  an 
ethylene-air  mixture  at  P=60  torr  and  0=0.5. 


ii.  Pure  Rotational  Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS)  studies  of  oxidation 
and  heat  release  in  ethylene-air  nanosecond  pulsed  discharges. 


While  measurement  of  the  temporal 
evolution  of  species  mole  fraction  after 
initiation  of  a  single  discharge  pulse,  sueh  as 
those  summarized  in  Figs.  2  and  3  above  of 
have  served  to  validate  modeling  predictions 
for  key  low  temperature  plasma  oxidation 
processes,  because  the  energy  coupled  to  the 
plasma  by  a  single  nanosecond  pulse  is  fairly 
low  (of  the  order  of  -0.1  mJ/cnr /pulse  or 
-0. 1  mcV/molcculc/pulsc),  detectable  heat 
release  requires  from  a  few  tens  to  a  few 
hundred  pulses  generated  at  a  high  pulse 
repetition  rate  (-10-100  pulses/msee).  To 
address  this  need  a  major  activity  of  the 
current  program  has  been  to  assemble  a  new, 
to  our  laboratory,  pure  rotational  Coherent 
Anti-Stokes  Raman  Spectroscopy  (CARS) 
instrument  and  to  use  it  for  a  scries  of 
measurements  of  rotational-translational 
temperature  rise  in  air,  ethylene-air,  and 


T,  °C 


Figure  4:  Comparison  of  experimental  spatially  averaged 
(emission)  and  centerline  (CARS)  temperatures  in  a 
repetitively  pulsed  nanosecond  discharge  in  air  at  P=40 
torr  with  plasma  ehemistrv  model  prediction. 
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hydrogen-air  mixtures,  excited  by  “bursts”  of  10  to  approximately  600  high  voltage,  nanosecond 
duration  pulses,  at  40  kHz  pulsed  repetition  rate. 

Figures  4-6  summarize  the  most  important  findings  for  ethylene-air  mixtures,  more  details  of 
which  can  be  found  in  publication  .  Figure  4  shows  time  resolved  air  plasma  temperatures 
inferred  from  both  N:  second  positive  emission  spectroscopy  and  CARS,  along  with  plasma 
chemistry  model  predictions  for  a  pulse  burst  discharge  at  P=40  torr  and  40  kHz  pulse  repetition 
rate.  Both  emission  (i.e.  spatially 
averaged)  and  CARS  (centerline) 
temperatures  are  very  close  to  each  other 
until  approximately  7.5  msec  (300  pulses). 

After  this  time,  the  emission  temperature 
rise  becomes  slower  and  starts  to  level  off, 
indicating  significant  heat  transfer  losses, 
while  the  CARS  temperature  continues  to 
increase  up  to  17  msec  (680  pulses).  At 
t=  1 5  msec  (600  pulses),  the  difference 
between  CARS  and  emission  temperature 
is  about  100°  C.  Note  that  the  plasma 
kinetic  model  assumes  a  constant  wall 
temperature  of  Tw=300  K,  which  appears 
to  adequately  describe  the  transverse 
temperature  distribution  in  the  air  plasma. 

Figure  5  compares  the  CARS  temperature  with  the  centerline  temperature  predicted  by  the 
kinetic  model  in  air,  lean  (cp=0. 1),  and  stoichiometric  (cp=1.0)  ethylene-air  mixtures  at  P=40  Torr 
and  a  burst  repetition  rate  of  40  kHz.  Both  experimental  and  predicted  temperatures  arc  plotted 
as  functions  of  discharge  burst  duration.  A  few  general  trends  are  clearly  observed.  First, 
comparison  of  experimental  temperatures  in  air  and  ethylene-air  show  that  the  rate  of 
temperature  rise  in  both  fuel-air  mixtures  is  much  faster  than  in  air  during  the  first  5  msec  (200 
pulses),  approximately  350  '  C  verses  125  C.  Second,  the  CARS  data  show  that  during  this 
stage  the  rates  of  heat  release  in  lean  and  stoichiometric  ethylene-air  mixtures  are  very  close  to 
each  other,  although  the  fuel  mole  fractions  differ  by  a  factor  of  10.  Third,  experimental 
temperatures  in  air  and  in  the  lean  ethylene-air  mixture  reach  a  quasi  steady-state  at  t~15  msec 
(600  pulses),  while  temperature  in  the  stoichiometric  ethylene-air  mixture  continues  to  rise. 

Two  of  these  trends  are  well  reproduced  by  the  kinetic  model.  In  particular,  the  model  predicts  a 
significantly  more  rapid  centerline  temperature  rise  in  air-fuel  mixtures  than  in  air,  consistent 
with  the  CARS  data.  Also,  predicted  rates  of  centerline  temperature  rise  in  air-fuel  mixtures  at 
<p=0. 1  and  1 .0  are  very  close  to  each  other  until  t=5  msec.  However,  the  model  predicts  a 
significantly  lower  temperature  in  the  stoichiometric  air-fuel  mixture  at  t>5  msec  than  is 
measured  experimentally.  Specifically,  the  model  predicts  that  by  5  msec  (200  pulses)  the 
temperature  reaches  a  quasi  steady-state,  whereas  the  CARS  data  indicates  that  it  continues  to 
rise.  To  explain  these  general  trends,  it  is  noted  that  the  rate  of  heat  release  in  the  plasma  as  a 
function  of  number  of  pulses  in  the  burst  depends  significantly  on  mixture  composition.  Figure 
6,  which  plots  atomic  oxygen  and  ethylene  mole  fractions  predicted  by  the  plasma  kinetic  model 
in  air,  and  (p=0.1  and  1.0  mixtures  (i.e.  at  the  conditions  of  Fig.  2),  provides  some  insight  into 
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Figure  5:  Comparison  of  experimental  centerline  (CARS) 
temperatures  in  a  repetitively  pulsed  nanosecond  discharge 
in  air  and  in  ethylene-air  (0=0. 1  and  O  =1 .0)  at  P=40  torr 
with  the  plasma  chemistry  model  prediction. 
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these  trends.  In  the  absence  of  fuel,  O  atoms 
formed  during,  and  immediately  after,  the 
discharge  pulse  decay  relatively  slowly  in 
recombination  reactions,  O  +  O2  +  M  — »  O3  + 

M  and  O  +  O3  — *  2O2.  In  air,  the  relatively 
low  rate  of  recombination  results  in  a 
considerable  build-up  of  O  atoms,  to  mole 
fractions  of  up  to  ~1%,  which  is  comparable 
to  the  0.7%  initial  mole  fraction  of  ethylene  in 
the  lean  air-fuel  mixture  (q>=0. 1 ).  However, 
in  ethylene-air  mixtures,  O  atoms  rapidly 
react  with  ethylene  due  to  the  parallel 
processes,  O  +  C2H4  — >  CHi  +HCO  and  O  + 

C2H4  — >  H  +  CH2CHO,  both  of  which  are 
relatively  fast  at  low  temperatures.  These 
low-temperature  reactions  of  atomic  oxygen 
with  ethylene  initiate  a  number  of  net  exothermic  hydrocarbon  reactions,  resulting  in  accelerated 
temperature  rise,  whieh  the  model  accurately  predicts,  as  can  be  seen  in  Fig.  5. 

Because  of  these  rapid  reactions  with  ethylene,  O  atom  number  densities  in  ethylene-air  mixtures 
are  much  lower  than  in  air  under  otherwise  identical  conditions.  This  can  be  seen  by  comparison 
of  the  solid  red  (air)  and  black  (ethylene-air,  0=1)  curves  in  Fig.  6.  Further,  as  long  as  fuel  is  in 
excess,  the  rate  of  chemieal  energy  release  (and  therefore  the  rate  of  temperature  rise)  in  air-fuel 
mixtures  is  dictated  by  the  rate  of  O  atom  generation,  whieh  is  only  slightly  affected  by  the 
ethylene  mole  fraction  in  the  mixture  (0.7-6. 5%).  For  this  reason,  the  predieted  initial  heating 
rates  in  the  lean  and  stoiehiometrie  ethylene-air  mixtures  are  very  close  to  each  other,  which  is 
again  consistent  with  the  experimental  CARS  temperature  measurements  shown  in  Fig.  4.  After 
nearly  all  fuel  in  the  lean  mixture  has  been  oxidized  (after  -200  discharge  pulses,  or  -5  msee), 
the  heating  rate  is  predicted  to  slow  down,  and,  simultaneously,  the  O  atom  number  density  is 
predicted  to  start  rising  rapidly.  However,  sinee  the  stoichiometric  mixture  has  ten  times  more 
fuel,  its  mole  fraction  remains  approximately  constant  (dashed  black  curve  in  Fig.  3),  as  does  the 
corresponding  predicted  heating  rate.  For  this  reason,  the  stoiehiometrie  mixture  continues  to 
generate  additional  heat  as  the  number  of  pulses  increases  beyond  -200-400  (5-10  msee).  The 
most  likely  reason  for  the  significant  difference  between  experimental  and  predieted  centerline 
temperatures  in  the  stoiehiometrie  mixture  after  large  number  of  pulses  (-400  pulses,  or  -10 
msee)  is  heating  of  the  test  section  walls.  At  these  conditions,  the  use  of  the  constant  wall 
temperature  boundary  condition  (Tw=300  K)  is  likely  no  longer  justified,  since  the  rate  of  heat 
transfer  from  the  plasma  to  the  eold  wall  is  significantly  overestimated. 

This  work  is  described  in  more  detail  in  Publications  5  and  14. 

iii.  Single  photon  Laser  Induced  Fluorescence  (L1F)  studies  of  NO  creation  and  lost  kinetics 
in  air,  methane/air,  and  ethylene/air  nanoseeond  pulsed  discharges. 

While  NO  production  and  destruction  in  non-equilibrium  plasmas  has  traditionally  been  studied 
from  the  perspective  of  pollutant  mitigation,  it  can  also  play  a  role  in  low  temperature  oxidation 
due  to  the  process  HO2  +NO  — *•  NO2  +  OH.  This  provides  both  a  low  temperature  means  of 
producing  OH,  with  subsequent  highly  exothermic  H  abstraction,  RH  +  OH  — »  R  +  H2O,  and  a 
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Figure  6:  Plasma  chemistry  model  predictions  for 
O  (solid  curves)  and  ethylene  (dashed  curves) 
mole  fractions  for  air  and  ethylene -air  mixtures. 


means  to  mitigate  a  key  low  temperature  ehain  termination  proeess.  Indeed,  it  has  been 
observed  in  equilibrium  low  temperature  flow  reactor  studies  by  Bromley,  et  al.,[3],  that  traee  (a 
few  ppm)  quantities  of  NO  aeeelerate  low  temperature(~650-700  K)  oxidation  of  butane. 

To  study  the  implications  of  this  for  nonequilibrium  PAC,  NO  studies  were  performed  using  the 
well  known  single  photon  X  (2II)  — >  A  (:X+)  LIF  transition  at  -  226  nm.  Note  that  this  is  very 
elose  to  the  atomie  oxygen  TAL1F  absorption,  enabling  the  use  of  an  essentially  identical 
experimental  apparatus.  Principal  results  of  these  studies  in  air  are  summarized  in  Figs.  7  and  8. 
Figure  7  shows  experimental  data  (black  circles),  as  well  as  modeling  predictions  for  NO  and  a 
number  of  other  key  species  using  a  “baseline’'  air  plasma  model  consisting  of  a  set  of  ordinary 
differential  equations  for  number  densities  of  a  large  number  of  neutral  and  charged  species 
produced  in  the  plasma,  as  well  as  the  energy  equation  for  predicting  the  temperature  of  the 
mixture.  This  set  of  equations  is  coupled  with  the  steady,  two-term  expansion  Boltzmann 
equation  for  the  electron  energy  distribution  function  (EEDF)  of  the  plasma  electrons,  using 
experimental  eross  sections  of  electron  impaet  electronic  excitation,  dissociation,  ionization,  and 
dissociative  attachment  processes.  Full  details  are  provided  in  publication  .  Focusing  on  the 
solid  blaek  curve  in  Fig.  7,  it  ean  be  seen  that  the  baseline  model  predicts  peak  NO  mole  fraction 
which  is  a  factor  of  -four  lower  than  that  observed,  and,  more  critically,  a  characteristic  time 
scale  for  peak  NO  concentration  which  is  longer,  by  -two  orders  of  magnitude,  than  that 
observed  experimentally  (-250  microseconds,  experimental,  vs  25  msec  from  baseline  model). 
Note  that  the  baseline  prediction  is  dominated  by  the  ordinary  equilibrium  Zeldovich  mechanism 

N2  +  0  — NO  +  N  (la) 

N  +  02-+N0  +  0  (lb) 


which  predicts  very  slow  NO  formation  at  low  temperature.  In  addition,  while  not  shown,  time- 
resolved  spontaneous  emission  measurements  showed  conclusively  that  the  N2  (C  II),  and 
N^A’X)  exeited  states  deeay  on  time  scales  of  10-20  nsee,  and  1-2  mieroseeonds,  respectively 
so  that  signficiant  NO  formation  via  processes  such  as  N2(A'I)  +  O  — »  NO  +  N  can  also  not 
explain  the  observations  in  Fig.  7. 
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Fig.  7:  Experimental  O  atom  and  NO 
concentrations  (symbols)  and  speeies 
concentration  predicted  by  the  baseline  model 
(lines)  after  a  single-pulse  discharge  in  air.  P 
60  torr. 


Fig.  8:  Experimental  O  atom  and  NO 
eoneentrationd  (symbols)  and  speeies 
concentration  predieted  by  Model  11  (“N:(vy' 
model,  lines)  after  a  single-pulse  discharge  in  air. 
P=  60  torr. 


6 


After  a  detailed  analysis,  documented  in  Uddi,  et  al,  2009  (Publication  3),  it  was  determined  that 
the  most  likely  process  responsible  for  the  observed  rate  of  low  temperature  NO  formation  is  the 
reaction  of  highly  vibrational ly  excited  ground  electronic  state  N2  with  ground  state  atomie 
oxygen,  ie 

N2(X,v>12)  +  0(3P)->N0  +  N(4S)  (2) 

where  N2(X,v>12)  is  formed  by  eollisional  quenehing  of  as  N2  (A  ’l) 

N2  (A  3I)  +  02 ->N2(X,v)  +  02(b  1  E)  (3) 

The  predictions  of  this  “N2  (X,v)”  model  are  shown  in  Fig.  8. 

Figures  9  and  10  summarize  results  in  tp=0.5  methane/air  (Fig  9)  and  ethylene/air  (Fig  10) 
mixtures,  where  experimental  data  is  compared  with  the  baseline  and  N2  (X,v)  models 
(proeess(2)),  as  well  as  a  third  model, 

N(4S)  +  02(b'l)  ->  NO  +  0(3P)  (4) 

where  02(b'Z)  is  formed  via  process  (3),  above.  It  can  be  seen  that  for  both  fuel/air  mixtures,  the 
02(b'S)  provides  better  agreement  with  the  experimental  data,  although  in  both  eases  the 
predictions  of  the  O^b'E)  and  N2  (X,v)  models  arc  quite  similar.  Note  also  that  the  experimental 
results  for  (p=0.5  methane/air  are  quite  similar  to  that  for  pure  air,  whereas  peak  NO  production 
in  (p=0.5  ethylene/air  is  lower  by  a  faetor  of-two. 


NO  mole  fraction 


Time,  seeonds 


Figure  9.  NO  concentrations  after  a  single-pulse 
discharge  in  methane-air  for  the  three  models 
used.  P=60  torr,  equivalence  ratio  (p~0.5,  pulse 
energy  1.0  mJ. 


NO  mole  fraction 


Time,  seeonds 

Figure  10.  NO  concentrations  after  a  single-pulse 
discharge  in  ethylene-air  for  the  three  models 
used.  P=60  torr,  equivalence  ratio  (p=0.5,  pulse 
energy  1.0  mJ. 


This  work  is  described  in  more  detail  in  Publications  4  and  12. 
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iv.  Pure  rotational  Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS)  studies  of  oxidation, 
heat  release,  and  ignition  in  hydrogen-air  nanosecond  pulsed  discharges. 

Furthering  our  understanding  of  plasma  assisted  hydrogen/air  combustion  is  clearly  important 
due  to  its  potential  utility  for  application  to  future  high  speed  (ie  SCRAMJET)  propulsion 
systems.  Its  inherent  simplicity  also  makes  it  an  ideal  platform  for  PAC  code  development  and 
validation.  During  the  final  year  of  the  three  year  program  of  research  summarized  in  this  report, 
a  significant  effort  has  gone  into  studies  of  hydrogen'air  PAC  including:  (i).  Experimental  CARS 
measurements  of  oxidation  and  heat  release,  (ii),  hydrogen/air  plasma-ehemieal  code 
development,  and  (iii),  hydrogen/air  PAC  sensitivity  analysis.  Each  of  these  elements  will  be 
summarized  below. 

Hvdrotzen/Air  Plasma-Chemical  Model  Development 

The  new  hydrogen/air  plasma-ehemieal  code  is  described  in  detail  in  Zuzeek,  et  al.,  2010  (Pub 
14)  and  will  be  briefly  summarized  here.  The  nanosecond  discharge  portion  of  the  model  is 
patterned  after  that  used  previously  for  air,  methane-air,  and  ethylene-air  discharges  (Topics  i-iii, 
above).  Briefly,  the  present  model  incorporates  a  nonequilibrium  air  plasma  chemistry  model 
developed  by  Kossy,  et  al  [4],  expanded  to  include  hydrogen  dissociation  processes  in  the 
plasma,  and  the  hydrogen-oxygen  chemistry  model  (22  reactions  among  H,  O,  OH,  Hj,  O2,  H2O. 
HO2,  and  H2O2)  developed  by  Popov  [5],  Species  concentration  equations  are  coupled  with  the 
two-term  expansion  Boltzmann  equation  for  the  energy  distribution  function  of  plasma  electrons, 
w  ith  electron  impact  cross  sections  taken  from  the  literature.  The  full  list  of  air  plasma  processes 
incorporated  into  the  model  and  their  rates  is  given  in  a  Publication  3.  The  model  also 
incorporates  results  from  a  new,  first  principles  model  describing  energy  coupling  to  the  plasma 
from  nanosecond  discharge  pulses,  which  will  be  described  separately  under  program  element  v, 
below,  and  quasi-one-dimensional  conduction  heat  transfer,  which  for  the  eases  to  be  discussed 
in  the  findings  section  assumes  a  constant  wall  temperature  boundary  condition  of  Tw  =  300  K. 
This  boundary  condition  can  be  modified  as  appropriate. 

Experimental  CARS  Results 

Figures  1 1  and  12  summarize  the  most  important  finding  in  hydrogen-air  mixtures.  Figure  1 1  is 
analogous  to  Fig.  5,  showing  experimental  CARS  temperatures  as  a  function  of  number  of  pulses 
in  a  40  kHz  burst  in  40  Torr  of  air,  and  hydrogen-air  at  a  variety  of  equivalence  ratios,  along  with 
modeling  predictions.  Unlike  Fig.  5,  however.  Fig.  1 1  shows  that  at  burst  durations  in  the 
approximate  range  16-20  msec  (640  -  800  pulses)  the  experimental  temperatures,  as  well  as  the 
model  predictions,  for  both  <p  =  0.5  and  tp  =1  exhibit  a  distinct  maximum,  and  then  decay. 
Furthermore,  as  shown  in  Fig.  12,  the  O2  mole  fractions  for  the  stoichiometric  mixture,  obtained 
from  the  Sandia  CARS  fitting  code  used  for  inference  of  rotational  temperature  from  the 
experimental  CARS  spectra  and  which  floats  this  parameter,  decreases  very  rapidly  in  the  region 
of  15  -  20  msec  duration  bursts.  Figure  12  also  shows  that  O2  mole  fraction  decay  agrees  well 
with  the  plasma  kinetic  model  predictions.  The  observed  temperature  maximum,  coupled  with 
the  observed  rapid  loss  in  O2,  is  a  strong  indication  that  ignition  has  occurred  in  these  mixtures. 
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Number  of  pulses 
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Figure  1 1:  Comparison  of  experimental  centerline  (CARS)  temperatures  in  a  repetitively  pulsed  nanosecond 
discharge  operated  in  repetitive  burst  mode  in  air  and  in  hydrogen-air  at  4>- 0.05,  0.5,  and  1.0  at  P=40  Torr  with 
plasma  chemistry  kinetic  model  predictions.  Pulse  repetition  rate  is  v=40  kHz,  burst  repetition  rate  10  Hz. 


Mole  fractions 


Time,  msec 


Figure  12:  Oxygen  mole  fraction  inferred  from  CARS  synthetic  spectra  for  hydrogen-air  at  the  conditions 
of  Fig.  4  (for  <)>=1.0).  Also  shown  are  hydrogen  and  oxygen  mole  fractions  predicted  by  a  plasma  ehenneal 
kinetic  model  for  the  same  conditions.  Oxygen  depletion  during  plasma  chemical  oxidation  /  ignition  is 
apparent. 
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Sensitivity  Analysis 


In  addition  to  hydrogcn/air  PAC  code  validation  summarized,  an  extensive  sensitivity  analysis 
has  been  performed,  which  is  described  in  detail  in  Zuzeek,  et  al  (2010)  and  Bowman,  et  al. 
(2010).  Summarizing  the  results  of  that  analysis,  very  briefly,  it  is  found  that  for  single 
discharge  pulse  oxidation,  initial  low  temperature  heating,  in  excess  of  that  resulting  from  direct 
energy  input  from  the  discharge,  is  controlled  primarily  by  the  following  sequence  of  three 
processes,  all  of  which  are  relatively  rapid  at  300  K, 


H  +  02  +  M 

ho2  +  m 

(5) 

o  +  ho2 

-»  OH  +  02 

(6) 

OH  +  H2 

h  +  h2o 

(7) 

Note  that  under  these  highly  non-equilibrium  conditions,  where  atomic  oxygen  is  formed  in 
much  higher  quantities  than  in  low  temperature  equilibrium  combustion,  the  three  body 
recombination  process  (5)  actually  acts  as  a  source  of  OH  and  subsequent  heat  release,  as 
opposed  to  a  radical  sink.  At  intermediate  temperatures,  500  -  600  K,  the  rate  of  fuel  oxidation 
increases  significantly,  primarily  due  to  four  additional  reactions  controlling  chain  branching. 


H  +  H02 

-*■  h2o  +  0 

(8) 

h  +  ho2 

-»  OH  +  OH 

(9) 

H  +  H02 

^h2  +  o2 

(10) 

O  (3P)  +  H2— >  H  +  OH. 

(11) 
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Figures  13  and  14  summarize  findings  for 
burst  mode  oxidation,  heating  and  ignition. 
Figure  13  compares  predicted  temperature 
rise  for  the  conditions  of  Fig.  11  (at  cp=  1 )  for 
the  full  plasma  kinetic  model  (solid  curve) 
and  a  reduced  chemistry  set  that  includes 
processes  (5)  -  (11)  above,  plus  the 
following  additional  two, 

H  +  02-»0  +  0H  (12) 

OH  +  HO:  -»  H20  +  02  (13) 

It  can  be  seen  that  the  predictions  are 
essential  identical.  Figure  14  is  analogous 
to  Fig.  13,  but  shows  the  temporal  evolution 
of  key  species  during  burst  mode 
oxidation/ignition  predicted  by  the  full  (left) 
and  reduced  (right)  chemical  model.  It  is 
again  found  that  the  two  sets  of  predictions 
arc  essentially  identical. 


▲ 


0  5  10  15  20  25  30 

Time,  msec 

Figure  13:  Temperature  during  Hi-air  oxidation/ignition  in 
burst  mode  at  v=40  kHz,  predicted  by  full  and  reduced 
reaction  sets.  Hi-air,  ()>=1,  P=40  torr.  initial  temperature 
T=300  K.  CARS  temperatures  shown  for  comparison. 
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Figure  14:  Species  mole  fractions  during  burst  mode  hydrogen-air  oxidatioa'ignition,  predicted  by  a  plasma 
chemical  kinetic  model  at  the  conditions  of  Fig.  13.  Full  reaction  set  (left)  and  reduced  set  (right). 


Note  that  the  results  summarized  above  are  presented  in  more  detail  in  Publieations  15  and  16.  A 
manuscript  has  also  been  submitted  to  the  2010  International  Symposium  on  Combustion,  to  be 
held  in  Beijing,  China  in  August  of  2010  (Publication  8). 

v.  Development  of  a  new,  non-empirieal  model  for  nanosecond  pulsed  discharge  plasma 
dynamics  and  energy  coupling. 

Nanosecond  pulsed  discharge  modeling  has  focused  on  two  closely  related  problems,  (a)  The 
temporal  dynamics  of  nanosecond  pulsed  discharge  development  and  (b)  Energy  coupling  to  the 
plasma  as  a  function  of  heavy  species  density,  which  vanes  with  discharge  pulse  number  when 
operating  in  burst  mode.  Both  analytic  and 
numeric  models  have  been  developed,  some 
principal  results  of  which  are  summarized  below. 

Full  details  are  provided  in  Publication  7. 

Figure  15  illustrates  the  one-dimensional 
geometry  that  is  used  for  the  discharge  model. 

Referring  to  Fig.  1,  the  discharge  section  can  be 
modeled  as  a  pair  of  electrodes,  separated  by 
distance  L  =  1  cm,  covered  with  a  dielectric  layer 
with  dielectric  permeability  e  (the  quartz  channel). 

As  illustrated  in  Fig.  16,  the  applied  voltage  is 
assumed  to  have  a  Gaussian  temporal  distribution 
with  full  width  at  half  maximum  of  -25  nsec  and 
peak  voltage  of  -18  kV.  Skipping  the  details,  the 
model  predicts  that  during  the  initial  stage  of  the 
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Figure  15.  Schematic  of  the  discharge  model  geometry 
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pulse,  when  the  applied  voltage  is  less  than  the  breakdown  voltage,  the  clcetric  field  between  the 
electrodes  is  uniform,  equal  to  Vapp/L.  When  the  applied  voltage  achieves  breakdown  (at  ~88 
nsec  in  Fig.  16),  a  high  field  sheath  (of  thickness  ls  ~  1  mm)  rapidly  forms  in  the  vicinity  of  the 
cathode.  Simultaneously  the  electric  field  within  the  “plasma”  drops  rapidly,  to  ~10%  of  the 
applied  voltage.  In  essence,  after  breakdown,  nearly  the  entire  applied  voltage  is  dropped  across 
the  sheath  region.  Figure  16  (left)  shows  that  with  the  exception  of  the  small  dip  in  the 
numerical  model  at  the  point  of  breakdown,  the  analytic  and  numeric  solutions  are  essentially 
identical. 


Plasma  electric  field,  kV/cm  Power  density,  kW/cm-' 


Figure  16.  Comparison  of  time-dependent  electric  field  (left)  and  coupled  power  (right)  in  the  plasma 
predicted  by  numerical  and  analytic  discharge  models,  along  with  applied  voltage  waveform. 


As  ean  be  seen  in  Fig.  16  (right),  the  majority  of  the  energy  coupled  to  the  plasma  occurs  in  a 
narrow  time  window  (~l-2  nanoseconds)  during  breakdown.  More  quantitatively,  the  model 
predicts  that  the  energy  coupled  to  the  plasma  during  the  pulse  is  limited  primarily  by  charge 
accumulation  on  the  dielectric  surfaces  covering  the  electrode  plates.  Basically,  the  coupled 
pulse  energy,  Qcap,  is 

q  »— c  V2  =i3eA  y2 

where  Ctotai  is  the  total  capacitance  of  the  charged  dielectric  plates,  e  is  the  dielectric 
permeability,  A  is  the  electrode  surface  area,  1  is  the  dielectric  plate  thickness,  and  Vo  is 
breakdown  voltage.  Note  that  since  breakdown  voltage  in  the  plasma  predicted  by  the  model, 
Vo=8-10  kV  at  the  conditions  of  the  present  experiments,  is  significantly  lower  than  the  peak 
pulse  voltage  applied  to  the  electrodes,  Vapp=20  kV  (Fig.  16  left),  the  coupled  pulse  energy  is 
only  ~  1 5-25%  of  the  maximum  energy  stored  in  the  capacitance  of  the  load  The  remainder  is 
reflected  back  to  the  power  supply  when  the  pulsed  voltage  is  reduced  to  zero.  It  ean  be  seen 
that  the  pulse  energy  coupled  to  the  plasma  can  be  significantly  increased  if  the  thickness  of  the 
dielectric  layers  is  reduced. 
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The  new  model  provides  a  theoretical  foundation  for  predicting  totaled  coupled  pulse  energy  as  a 
function  of  pulse  number  in  a  pulse  burst.  It  is  shown  in  publication  6  that  the  totaled  coupled 
energy  per  molecule  is  approximately  constant,  so  that 

1  oo  V2  1 

q„.  » - — - «  0.28  meV  /  molecule  (13) 

pu,se  2  21 L  Ne 

In  other  words,  as  the  temperature  increases  during  a  burst  (at  constant  pressure),  the  number 
density  of  the  plasma  drops,  resulting  in  a  linear  drop  of  total  coupled  pulse  energy  with 
temperature.  Physically  this  occurs  because  the  drop  in  density  results  in  a  drop  in  breakdown 
voltage,  reducing  the  total  coupled  pulse  energy  according  to  (12). 
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